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Microwave probes Dipole Blockade and van der Waals Forces in a Cold Rydberg Gas 
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We show that microwave spectroscopy of a dense Rydberg gas trapped on a superconducting atom 
chip in the dipole blockade regime reveals directly the dipole-dipole many-body interaction energy 
spectrum. We use this method to investigate the expansion of the Rydberg cloud under the effect 
of repulsive van der Waals forces and the breakdown of the frozen gas approximation. This study 
opens a promising route for quantum simulation of many-body systems and quantum information 
transport in chains of strongly interacting Rydberg atoms. 

PACS numbers: 03.67.-a,32.80.Ee,32.30.-r 


The strong dipole-diplole interaction between cold Ry¬ 
dberg atoms is the focus of an intense theoretical and ex¬ 
perimental interest [1-3] . It provides an efficient platform 
for quantum information processing, with quantum gates 
based on the manipulation of single trapped atoms [4, 5] . 
It leads to optical non-linearities [6, 7], even at the single 
photon level [8-11]. Finally, it opens the way to quan¬ 
tum simulators of self-organization processes and phase 
transitions [12-16]. 

A variety of methods have been implemented to in¬ 
vestigate the dipole-dipole interaction. An early ob¬ 
servation of a Rydberg excitation line broadening at 
high density provided a direct evidence of the interac¬ 
tion [17]. Dipole blockade [18-21] was demonstrated 
with individually trapped atoms [22], observed on Ry¬ 
dberg atom counting statistics in atomic ensembles [23- 
26] or on collective single-photon emission [27]. Spatial 
self-organization has been evidenced by individual atom 
imaging in gases [28] and in optical lattices [29]. 

We present here a direct measurement, based on mi¬ 
crowave spectroscopy, of the interaction energy distribu¬ 
tion in a cloud of ultra-cold Rydberg atoms in the block¬ 
ade regime. We use it to observe the atomic cloud ex¬ 
pansion driven by the repulsive van der Waals (vdW) 
forces, revealing the limits of the frozen gas approxima¬ 
tion [30]. The observations are in good agreement with 
Monte Carlo simulations. 

The core of the set-up is sketched on figure 1(a) (de¬ 
tails in [31, 32]). Rubidium atoms are trapped on a su¬ 
perconducting atom chip, which is cooled down to 4.2 K. 
In the experimental sequence, they first effuse from a 
2D-Magneto Optical Trap (MOT) placed at room tem¬ 
perature outside the cryostat. They are then caught a 
few millimeters away from the reflecting front chip sur¬ 
face in a mirror-MOT at a temperature T « 100 /rK. 
The MOT quadrupolar magnetic held is generated by 
centimeter-sized superconducting coils. The atoms are 
then transferred into a compressed mirror-MOT located 
700 ^m away from the chip. Its magnetic held is the 
superposition of that produced by the [/-shaped super- 
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FIG. 1. (color online) (a) Scheme of the central part of the 
setup, (b) Scheme of the superconducting atom chip. The 
axes are dehned in the two panels. The origin O is at the 
center of the Z wire LG. 


conducting wire connecting pads J and L [Fig. 1(b)] with 
a uniform bias. 

The magnetic held is then transiently switched off 
while the atoms are cooled down to 12 /rK by an op¬ 
tical molasses. They are optically pumped into the 
5)S,F = 2,tof = 2 level and transferred into a lohe- 
Pritchard trap. Its held is a superposition of that gen¬ 
erated by the Z-shaped GL wire (width: 70 ^m) with 
a uniform bias. The held at the bottom of the trap is 
aligned along the x quantization axis (dehned in Fig. 1). 
The atoms are then evaporatively cooled just above the 
BEC transition using a radio-frequency held radiated by 
the KM wire [31]. The total cloud preparation time is 
about 8 seconds. The hnal thermal cloud, 300 p,m away 
from the chip surface, contains « 12 000 atoms at a tem¬ 
perature of about 500 nK. The Gaussian atomic density 
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peaks at 2.4 10^^ cm“^, with dispersions 2cra; = 44 
and 2ay ~ 2az = 11 /im along the three axes. This mod¬ 
erate density makes the interaction between the Rydberg 
atoms and the ground state cloud (a few tens of kHz [33]) 
negligible as compared to the vdW interaction between 
Rydberg atoms. 

A large fraction of the atoms released from the trap is 
adsorbed on the gold front chip surface. The deposited 
atoms create large stray electric field gradients in the 
trap. This detrimental effect is avoided by coating the 
full chip surface with a metallic rubidium layer, stable 
at 4 K, which increases dramatically the Rydberg levels 
coherence time [32]. 

We induce the two-photon 5S —>■ bos' transition by 

2 fis 780 nm ‘red’ and 480 nm ‘blue’ laser pulses. The 
detuning with respect to the intermediate 5 P 3/2 level is 
540 MHz. The red and blue lasers propagate along the x 
axis [Fig. 1(a)], with 150 and 22 ^m waists, 50 fiW and 
8 mW powers, and a+ and a- polarizations respectively. 
They excite the 60S'i/2, mj = 1/2 sublevel with a 400 kHz 
two-photon Rabi frequency. Ten excitation pulses, at a 

3 ms time interval, are sent on the atomic cloud without 
noticeable heating. Both lasers are frequency-locked to a 
transfer cavity, whose length is stabilized to a rubidium 
saturated absorption line. At low Rydberg density, the 
excitation linewidth is 7 = 600 kHz, after minimization 
of the residual electric field in the y direction using elec¬ 
trode / facing the chip (kept at 0 V) [Fig. 1(a)]. The 
contribution of the residual electric field gradients is be¬ 
low 50 kHz [32]. 

We detect the Rydberg atoms by field-ionization [Fig. 
1(a)], with a 90±10% detection efficiency [32]. An elec¬ 
tric field ramp, produced by electrode I, reaches at dif¬ 
ferent times the ionization thresholds of the 60S and 
57S levels involved in the spectroscopy (37 V/cm and 
41 V/cm respectively). The resulting ions are acceler¬ 
ated and deflected (electrode D) towards a channeltron 
counter (Sjuts Optotechnik KBL lORS-EDR). For short 
delays between laser excitation and detection, the ion sig¬ 
nal is broadened, an effect we attribute to state mixing 
in the field ramp in the presence of dipole-dipole inter¬ 
actions. We thus apply the ionizing field 150/iS after the 
laser pulse, leaving time for the Rydberg cloud to expand 
due to the repulsive vdW forces (see below). Atomic in¬ 
teractions are thus negligible during field ionization. A 
residual partial overlap of the ionization signals of 60S 
and 57S results in a transfer rate offset, which is mea¬ 
sured and subtracted from the data. 

For interatomic distances, i?, larger than ~ 3 /im, 
the dipole-dipole energy shift for two 60S atoms reads 
hCgofio/. A numerical diagonalization of the Hamilto¬ 
nian provides a repulsive vdW interaction with Ceofio = 
137.5 GHz/im®. Note that Penning ionization [34, 35] is 
then negligible (the free ions collected at the onset of the 
field ionization ramp are ~ 1% of the Rydberg signal). 

The spatial distribution of the Rydberg atoms in the 



FIG. 2. (color online). Microwave spectra of the 60S — 57S 
two-photon transition for (a) A = 0 (b) A = 1 MHz, (c) A = 
2 MHz. The dots are experimental with statistical standard 
deviation error bars and the solid line results from a simple 
model. The origin of the frequency axis (thin vertical line) 
corresponds to the resonance position in a low-density cloud. 

sample and thus their interaction energy can be con¬ 
trolled by adjusting the detuning A of the blue laser with 
respect to the resonance frequency of non-interacting 
atoms [26, 34, 36, 37]. In a simple incoherent excitation 
model at resonance (A = 0), dipole blockade precludes 
the excitation of more than one Rydberg atom inside a 
sphere with radius Rt = (Ceo.eo/b)^'^® — 8 /rm. For a 
blue detuning (A > 0), after the non-resonant excitation 
of a first ‘seed’ Rydberg atom, the excitation of a second 
one is resonantly ‘facilitated’ on a spherical surface of ra¬ 
dius Rf = (Cgo.eo/^)^^® [36]. As additional atoms get 
excited, the facilitation surface increases. An avalanche 
process forms a cluster around the seed atom. If A > 7, 
Rf < Rb and the mean interatomic distance is reduced. 
The interaction energy per atom thus increases with A. 

We probe the interaction energy distribution by mi¬ 
crowave spectroscopy. We choose a 60S—l^nS two-photon 
transition, with a narrow (few kHz width) spectral line 
at frequency vq /2 for non-interacting atoms [32] . Atomic 
interactions shift and broaden this line. Let us consider 
first a pair of atoms at a distance R. A weak microwave 
couples the initial [bos', AOS') state only to [bos', nS”) and 
Ins', bos') states. The latter are resonantly coupled by 
the vdW interaction. The interaction matrix involves 
diagonal, hCeo.n/R®, and off-diagonal exchange terms, 
hAgQ n/R^- The eigenstates are symmetric and antisym¬ 
metric combinations of the bare states, with energies sep¬ 
arated by 2hAeo,n/R^- We thus generally get two ex¬ 
citation lines for a single pair. The situation is much 
simpler if we choose n = 57 {vq = 116.457 GHz). The 
exchange term, Agg 57 = —0.33 GHz /rm®, is then much 
smaller than Cgg^gr = —43.3 GHz/xm®. We get a sin¬ 
gle excitation line, at frequency [i/g -I- A:/(i?)]/2, with 
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Av{R) — (Ceo ,60 ~ C'eo.s?)/^^ — ^ where 

77 ~ 1.316. 

For N interacting atoms, an immediate generalization 
leads to a frequency shift for atom number i 

(assumed to be far away from all other atoms in state 
57S), with 

= (Ceo.eo — C'eo.s?) ^ 

l<j<N 

i^3 

= r]AE^'i/h, ( 1 ) 

where AE^'^'> is the vdW energy shift of atom i in the GO.S' 
state due to its interaction with all other Rydberg levels. 

We first probe the spectrum with a 1 p,s microwave 
pulse at an adjustable frequency, (t/q + Av)f2, applied 
immediately after laser excitation, so that atomic motion 
plays no role. We use three laser detunings, A = 0, 1 
and 2 MHz, resulting in about 80, 60 and 40 detected 
Rydberg atoms for each laser pulse respectively. Less 
than « 3 atoms are transferred into the 57S state. Di¬ 
rect interactions between these 57S atoms can be ne¬ 
glected. Since the vdW 57S-60S interaction is attractive, 
some 57S atoms (less than 50%) undergo Penning ioniza¬ 
tion [34]. This process, taking place after the microwave 
pulse, reduces the measured transfer rate but does not 
alter the spectrum shape. 

Figure 2 presents, for the three A values, the fraction 
of atoms detected in 575' as a function of the scaled mi¬ 
crowave frequency, Av/rj. The points are experimental 
and the solid lines result from a Monte Carlo simula¬ 
tion assuming successive incoherent Rydberg excitations 
from the ground state. At each iteration step, we se¬ 
lect randomly a ground state atom according to the trap 
geometry. We randomly choose whether this atom is ex¬ 
cited or not, taking into account the vdW energy shift 
due to previously excited atoms, the full laser linewidth 
and intensity profile. The number of iterations is chosen 
to reproduce the observed average Rydberg atom num¬ 
ber. The final microwave spectrum is a convolution of 
the interaction energy distribution with the microwave 
pulse linewidth. We average 50 to 200 Monte Carlo real¬ 
izations and adjust the vertical scale so that the area of 
the spectrum fits the experimental one. This procedure 
provides us with the solid lines in Fig. 2. They are in 
fair agreement with the experiment. 

The A = 0 spectrum maximum [Fig. 2(a)] is 
frequency-shifted by about the laser linewidth 7 . Hence, 
the interatomic distances are close to the dipole blockade 
radius R^. The high frequency tail is due to atoms in the 
cloud bulk, with several neighbors at a distance close to 
Rb- We infer from the simulation that the size of the Ryd¬ 
berg atom ensemble is three times larger than that of the 
ground state cloud. This large broadening results from 
dipole blockade. Since atoms cannot be excited at short 
distances, laser excitation favors the tail of the thermal 
cloud Gaussian distribution. This precludes the observa- 



FIG. 3. (color online) Rydberg atoms interaction energy spec¬ 
trum as a function of the delay r between atomic preparation 
and microwave probe (0.5, 10, 15, 20, 30, 50, 100 ps from 
bottom to top). The points are experimental (error bars not 
shown for clarity) and the lines result from the model (see 
text). The different spectra are shifted vertically by a fixed 
offset (0.15) for the sake of clarity. 


tion of a clear Rydberg atom number saturation with a 
thermal cloud. 

For A = 1 and 2 MHz [Fig. 2(b) and (c)], the to¬ 
tal Rydberg atom number is reduced due to the non¬ 
resonant excitation. Nevertheless, the interaction energy 
increases with A, corresponding to an increasing peak 
Rydberg density. The average shift is of the order of 2A, 
as expected from an energy conservation argument stat¬ 
ing that — N.hA. The agreement with 

the Monte Carlo simulation is good for A = 1 MHz. For 
A = 2 MHz, the simulation predicts a bimodal structure 
with a narrow component centered on A. This peak is 
due to atoms excited on the outer ‘facilitation’ surface, 
with an initial interaction energy hA unmodified since 
there is no further excitation in their vicinity. Our sim¬ 
ple model, which does not take into account coherent pair 
excitation processes, overestimates this contribution. 

We then investigate the Rydberg gas expansion due to 
the repulsive vdW forces between 60S atoms. Figure 3 
presents (points) measurements of the interaction energy 
spectra at different delays between the laser and the mi¬ 
crowave pulses for A = 1 MHz. For this detuning, the in¬ 
ternal potential energy per atom is of the order of 50 /iK, 
two orders of magnitude above the initial temperature. 
As T increases, the spectrum gets narrower and converges 
to that for a dilute cloud with negligible interactions. The 
theoretical lines in Fig. 3 result from the Monte Carlo 
excitation model, followed by a direct integration of the 
Newtonian equations of motion. The calculation takes 
into account the finite lifetime (210/rs [32]) of the 60S 
state, even though the atomic decay contribution to the 
interaction energy reduction is small (at most 30%). The 
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agreement with the experimental data is good, reinforc¬ 
ing our confidence in the excitation model. 



FIG. 4. (color online) Variation of the average interaction 
energy as a function of time for (a) A = 1 MHz and (b) 
A = 2 MHz. The points (connected by a solid line for visual 
convenience) are experimental. The thick dot-dashed line re¬ 
sults from our model. The dashed line corresponds to a single 
Rydberg atom pair expansion model. 

We plot in figure 4 the time evolution of the average 
interaction energy for A = 1 and 2 MHz (points and 
connecting solid lines). We observe, as expected, a decay 
of the vdW potential energy at the benefit of atomic ki¬ 
netic energy, faster when the initial energy is larger. The 
dot-dashed lines present the results of the Monte Carlo 
model, in fair agreement with the data. Since the model 
slightly underestimates the initial energy, particularly for 
A = 2 MHz, it predicts a slightly slower expansion. The 
experimental expansion is much slower than that of a 
single Rydberg atom pair with the same average initial 
energy (dashed lines in figure 4). The expansion pro¬ 
ceeds in an interesting hydrodynamic regime, requiring 
more detailed studies. 

We use microwave spectroscopy for a direct measure¬ 
ment of the vdW interaction energy distribution in a cold 
atom sample containing up to 80 interacting Rydberg 
atoms. The spectra provide direct evidence of the dipole 
blockade regime. Time-resolved measurements give ac¬ 
cess to the dynamical evolution of the Rydberg atom 
cloud, whose fast expansion sets limits to the frozen gas 
approximation. 

This new diagnostic method leads to a direct observa¬ 
tion of self-organization and dynamical phase transitions 
in Rydberg atom-based quantum simulators [15]. The 
strong confinement possible with a superconducting atom 


chip opens the way to the realization of unidimensional 
atomic clouds, with long Rydberg coherence times [32] 
and, hence, to quantum simulations of interacting spin 
chains [16, 38-40]. In particular, the strong exchange 
vdW interaction between 60S and 61S levels could be 
used to simulate quantum transport processes [41]. 
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